Valli et al. Extreme Physiology & Medicine 2013, 2:7 
http://www.extremephysiolmed.eom/content/2/1/7 



Extreme Physiology 
& Medicine 



RESEARCH Open Access 



Minute ventilation and heart rate relationship for 
estimation of the ventilatory compensation point 
at high altitude: a pilot study 

Gabriele Valli 1 , Mattia Internullo 1 , Alessandro M Ferrazza 1 , Paolo Onorati 1 , Annalisa Cogo 2 and Paolo Palange 1,3 



Abstract 

Background: The ventilatory compensation point (VCP) is an exercise threshold which has been used in the design 
of training programs in sports medicine and rehabilitation. We recently demonstrated that changes in the slope of 
the minute ventilation to heart rate relationship (Ai/f/AHR) can be utilized for estimation of the VCP during 
incremental exercise at sea level (SL). We hypothesized that in hypoxic conditions, such as high altitude (HA), 
VCP can be also reliably estimated by Al/f/AHR. 

Methods: At SL and on immediate ascent to HA (5,050 m), six healthy subjects (42 ± 14 SD years) performed a 
maximal incremental exercise test on a cycle ergometer; 0 2 uptake {Vo 2 ), C0 2 output {Vco 2 )i Ve, and HR were 
measured breath-by-breath. The Al/f/AHR method for VCP estimation was compared to the standard method 
using the ventilatory equivalent for C0 2 {Ve/V C o 2 ) and end-tidal PC0 2 (P ET C0 2 ). The Al/f/AHR slope values below 
(S-|) and above (S 2 ) VCP were computed by linear regression analysis. 

Results: A significant difference between S, and S 2 was observed, at SL and HA, for both the Al/ f /AHR and 
Ve/Vco 2 methods for VCP estimation. A good agreement between the two methods (Al/ f /AHR vs. V E /Vco 2 ) was 
found for both environmental conditions; the mean difference + 2 SD of Vo 2 at VCP (VCP-Vo 2 ) was -22 + 112 ml/ 
min at SL and 39 ± 81 ml/min at HA. The VCP-Vo 2 was significantly lower at HA compared to SL; in addition, S ] and 
S 2 mean values were significantly higher at HA compared to SL. 

Conclusion: At HA, VCP may be reliably estimated by the Al/ f /AHR method. 
Keywords: High altitude, Ventilatory compensation point, Exercise, Ventilation, Heart rate 



Background 

In normal individuals, minute ventilation ( Ve ) during 
exercise displays a more marked rapid increase when 
work rate exceeds the heavy-intensity domains, as a 
compensatory effect of the metabolic acidosis [1,2], and 
this threshold is called the ventilatory compensation 
point (VCP) [3]. The VCP has been utilized for the de- 
sign of training programs in sports medicine and re- 
habilitation [4,5]. Practical and, in some instances, 
simple approaches for estimation of the VCP have been 
evaluated in the past by different authors. Several 
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strategies have been proposed, based either on Ve [6] or 
heart rate (HR) responses [7,8]. 

Since HR is usually reasonably linear with respect to 
oxygen uptake (Vo 2 ) over the entire work rate range [9], 
we have argued that the changes in the slope of incre- 
ment in Ve over HR response (AVf/AHR) can be uti- 
lized to estimate VCP during incremental exercise at sea 
level (SL) [10,11]. The AV £ /AHR method has shown a 
good agreement with the standard method for VCP esti- 
mation (i.e., increased ventilatory equivalent for CO2 
Ve/V C q 2 ) and decreased end-tidal PC0 2 (PetC0 2 ), in- 
dependently of exercise modality [11] and protocol [10]. 
The AVf/AHR method therefore may be a potentially 
promising method to estimate VCP during field tests or 
in difficult environmental conditions. 
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At high altitude (HA), compared to SL, V E is greatly 
increased at any given work rate [12]. This may have a 
detrimental effect on exercise tolerance because it may 
lead to a significant reduction in ventilatory reserve [13]; 
however, reduced regional blood flow and oxygen diffu- 
sion limitation in the lungs and exercising muscles are 
also believed to play an important role in limiting 
exercise tolerance at HA [14]. The HR response at any 
given work rate is also increased at HA [15], although 
maximal HR appears to be reduced [13-16]. No data, 
however, are available on the slopes of the Al/e/AHR re- 
sponse to exercise at HA. 

In the present study, we wanted to establish the reli- 
ability of the AV^/AHR method for the estimation of 
VCP at HA, an environmental condition associated with 
an augmented ventilatory and cardiac response to exer- 
cise. We hypothesized that the AV^/AHR method at 
HA would be as reliable as the standard method for 
VCP estimation. If so, this could lead to the develop- 
ment of less expensive exercise equipment, i.e., without 
gas analyzers, that could readily be easily used in sports 
medicine and in high altitude studies where calibaration 
of gas analyzers is more difficult. 

Methods 

Subjects 

Six healthy non-smoking SL resident subjects (two fe- 
males, four males; Table 1) provided written informed 
consent to participate in the study. The procedures and 
protocols were approved by the National Research 
Council (CNR, Italy) as part of the Ev-K2-CNR research 
program and were conducted in accordance with the 
Declaration of Helsinki. 

Protocols 

The SL phase was conducted in Rome, Italy, at approxi- 
mately 60 m of altitude. On day 1, the subjects 
performed standard pulmonary function tests and got 
familiarized with the maximal incremental exercise test. 



Table 1 Subjects' baseline characteristics 



Baseline characteristic Value 



Age (years) 




41 ± 15 




Height (cm) 




1.74 ± 0.10 




BMI (kg/m 2 ) 




23.6 ± 2.9 




FVC at SL (1, 


%pred) 


4.70 ± 0.90, 103 ± 




FEV, at SL (1, 


%pred) 


3.77 ± 0.75, 109 ± 


16 


FVC at HA (], 


%pred) 


4.59 ± 0.99, 1 00 ± 


14 


FEV, at HA (1 


, %pred) 


3.79 ± 0.98, 1 06 ± 


17 



FVC, forced vital capacity; FEV, forced expiratory volume; %pred, 
percent predicted. 



On the subsequent day, the formal maximal incremental 
exercise test was conducted until exhaustion. 

The HA phase of the study was conducted at the CNR 
Pyramid Laboratory, Lobuche, Khumbu, Nepal (5,050 m, 
barometric pressure approximately 410 Torr). After trav- 
elling by plane through Kathmandu (1,340 m) up to 
Lukla (2,860 m), all subjects performed a 10-day trek to 
the Pyramid Laboratory with the following intermediate 
stops for acclimatization: 3 days at Namche Bazaar 
(3,450 m) and then 2 days at Pheriche (4,252 m). Symp- 
toms of acute mountain sickness and resting arterial O2 
saturation (Sp0 2 ) were monitored. After 1 day of rest at 
the Pyramid Laboratory, the experimental phase started. 
On day 2, the subjects performed standard pulmonary 
function tests. Calibration of the spirometer (photoelec- 
tric digital turbine, diameter 28 mm, resolution 4 ml, 
Micro Kit, COSMED, Rome, Italy) was performed prior 
to each test, using a 3-1 syringe. 

Incremental exercise test 

On day 3, each subject performed a maximal incremental 
exercise test on a cycle ergometer. The pedaling frequency 
was set at 60 ± 5 rpm by subjects following an electronic 
audio signal generated each second. The test protocol 
consisted (a) 2 min of rest, (b) 4 min of exercise at 20 W, 
(c) the incremental phase with work rate increments of 15 
W/min at HA and 25 W/min at SL, (d) 6 min of recovery 
pedalling at 20 W. 0 2 uptake ( Vo 2 , standard temperature 
and pressure, dry (STPD)), C0 2 output (Vco 2 > STPD), mi- 
nute ventilation (Ve, body temperature and pressure satu- 
rated BTPS), and end-tidal partial pressures for 0 2 and 
C0 2 (Pet0 2 , P ET C0 2 ) were obtained breath-by-breath 
and edited to exclude occasional outlying breaths (>±4 SD 
of the local mean) as a result of coughs, swallows, sighing, 
or gasping. A 10-s average was utilized for subsequent 
analysis. 

For VCP estimation by analysis of the AV^/AHR rela- 
tionship (VCP^ £ / HR ), we utilized a least squares regression 
analysis. Using a dedicated software, we applied a 'best fit' 
line to the data which extended from end-exercise back to 
the sub-maximal point at which the linearity was lost (S 2 
region); this breakpoint was compared with the standard 
approach for VCP estimation (VCPy £ yy co ) based on the 

Ve/Vco 2 relationship, where Ve started to change out of 
proportion of Vqo 2 an d PetC0 2 started to fall [17]. From 
the aforementioned breakpoint, a second best fit line (Si) 
was applied to the data extending from the end of the 
warm-up phase upwards into the exercise data. If a second 
breakpoint was discernible, the Si line was then divided in 
two different regions (Si" from the breakpoint to 
VCPy E / HR and Si' from the end of warm-up up to the 
observed breakpoint), and the HR value at this breakpoint 
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was then compared with the value observed at lactate 
threshold (6l, estimated by the 'V-slope method' [3] and 
supported by standard ventilatory equivalent and end- 
tidal gas tension criteria [17]). 

Equipment 

The incremental exercise test was performed on a mech- 
anically braked cycle ergometer (828E, Monark Exercise 
AB, Varberg, Sweden). Prior to each test, the zero-load 
setting on the ergometer was checked, and a calibration 
was performed with a 4-kg weight. Ventilatory and pul- 
monary gas exchange variables were measured breath-by 
-breath in all tests using a portable system (K4b 2 , 
COSMED) which has been previously validated at HA 
[18]; the accuracy of the telemetric system has been pre- 
viously established [18,19]. The system comprised a face 
mask, analyzer unit (containing O2 and CO2 gas ana- 
lyzers), heart rate monitor, and battery. The analyzer 
unit with battery pack, face mask, and tubing (weight 0.8 
kg) was attached to the subject with a harness and 
connected to a personal computer by an Ethernet cable 
connection. The face mask contained a turbine for 
measurement of volume and flow; calibration was 
performed with a 3-1 syringe (Hans Rudolph, Kansas 
City, MO, USA) over a range of different flow profiles. 
Respired gas, sampled continuously from a port within 
the turbine via a Nation polymer capillary (PermaPure©, 
Toms River, NJ, USA), was analyzed at 100 Hz using 
rapid-response O2 (polarographic) and CO2 (infrared) 
analyzers (mean response time 120 ms) which were 
automatically thermostated and compensated for ambi- 
ent variations in barometric pressure, humidity, and en- 
vironmental temperature. Analyzers' calibration was 
performed using two precision-analyzed gas mixtures 
spanning the respired range. The volume and gas con- 
centration signals were sampled and digitized every 10 
ms, and time-aligned, i.e., correcting for the transport 
delay between the turbine and gas analyzers and for the 
analyzer rise time [20]. HR was measured from a chest 
strip and recorded every breath. Sp0 2 was monitored 
non-invasively by finger pulse oximetry (Masimo Rad-5, 
Masimo Corporation, Irvine, CA, USA). 

Statistical analyses 

Differences among measured responses were determined 
by a Student's paired t test. Pearson's product-moment 
correlation coefficient (R ) was used to identify correla- 
tions between criterion variables. The level of statistical 
significance was set at P < 0.05. Group data are 
presented as mean ± SD. The limits of agreement be- 
tween the Ve/Vco 2 an d tne AV^/AHR methods for 
VCP estimation were evaluated by the Bland-Altman 
analysis [21], where the individual differences are plotted 



against their respective means. We proceeded with such 
a type of analysis if a significant linear correlation 
between methods was previously observed. The same 
statistical approach was also performed to compare the 
correspondence between Gl and the breakpoint between 
Si" and Si' in the AVe/AHR relationship, when the lat- 
ter was detectable. 

Results 

The main results of the incremental tests are illustrated 
in Table 2. At HA, compared to SL, a significant reduc- 
tion in V^peaks . WRpeak, HR peak , Sp0 2 pe a k, and 8 L was 
observed; on the contrary, V^peak at HA was appreciably 
higher. 

Figures 1 and 2 show AV E /A 0l , Ve/Vqo 2 > P E t0 2 , 
P ET C0 2 vs. Vco 2 i and AV^/AHR relationship in two 
representative subjects at SL (upper panels) and HA 
(lower panels). As shown in Figure 3, in all subjects, a 
breakpoint in AVe/AHR, which occurred at the VCP es- 
timated by the Ve / Vco 2 method, was clearly discernible 
both at SL and HA. No significant differences were 
found in V 0l measured at VCP ( V COl -VCP; Table 3) 
between methods utilized to identify the threshold 

(VCP-AVVAHR vs. VCP-Ve/V C02 )- 
The Bland-Altman analysis confirmed the agreement 

between the Ve/Vqo 2 and AVe/AHR methods in esti- 
mating VCP, both at SL and at HA (Figure 4, Table 3); 
the mean bias ± 95% confidence interval of the between- 
method differences were -22 ± 112 ml/min of V o 2 for 
the SL protocol and 39 ± 81 ml/min of V 02 for the HA 
protocol (Figure 4). 

Although in some instances a lower breakpoint in the 
Si region was observed (four out of six subjects at SL 
and one out of six at HA), no significant agreement 
between such a breakpoint and the 9 L was obtained in 
terms of HR and V o 2 values. Votryii -VCP was signifi- 
cantly lower at HA (approximately 500 ml of V o 2 less) 



Table 2 Results of the incremental exercise tests at HA 
and SL 





HA 


SL 


P value 


Sp0 2peak (%) 


72.8 ± 7.3 


98 ± 1 .4 


<0.001 


WR peak (W) 


138 ± 26 


208 ± 44 


<0.01 


^o 2pEak (ml/min) 


1,690 ±533 


2,227 ± 511 


<0.01 


V'aW (ml/min) 


2,060 ±570 


2,820 ± 780 


<0.01 


14 P eak (l/min) 


123.4 ± 33.3 


91.3 ± 28.6 


<0.01 


9l (ml/min) 


1,015 ± 222 


1,302 ± 241 


<0.01 


HRpeak (bpm/min) 


147 ± 13 


1 68 ± 1 2 


<0.001 
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Figure 1 Representative subject #6 (male, 56 years). Exercise response variables in a representative subject at SL {empty circle) and at HA 
(filled circle). The block arrow indicates VCP; the white arrow indicates the 6 L . A single breakpoint in the \4/HR at VCP was clearly identifiable, both 
at SL and at HA (see text for further comments). 
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Figure 2 Representative subject #1 (male, 36 years). Exercise response variables in a representative subject at SL {empty circle) and at HA 
{filled circle). The black arrow indicates VCP; the white arrow indicates the 6 L . In this subject, a second and less evident breakpoint l/f/HR that 
occurred at the HR value observed at 0 L was identifiable (see text for further comments). 
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Figure 3 Exercise l/ E /HR relationship in all subjects. Exercise V E /HR response in the six subjects studied at SL {circle) and at HA {square)at free 
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Table 3 VCP, S v and S 2 values at HA and SL 

HA SL P value 

VCP - Al/ f /AHR (ml/min) 1 ,325.3 ± 214.7 1 ,822.0 ± 371.6 <0.0S 

VCP - Ve/Vco 2 (ml/min) 1,315.1 ± 180.8 1,839.1 ± 359.5 <0.05 



S H ,m " /b P) 1.09 ± 0.33 0.76 ±0.32 <0.01 

S \/HR ( ' /b P ) 4 ' 8 ± 25 1 4 ± 0 - 6 <0 - 05 

Su „, 0.043 ± 0.005 0.025 ± 0.001 <0.01 

S 2 . ,. 0.1 14 ± 0.052 0.047 ±0.01 <0.05 

S 2 /S Ve/hr 4.5 ±1.9 2.0 ±0.8 NS 

S 2 /S, 2.7 ±1.5 1.9 ±0.4 NS 



NS, not significant. 

independently of the method utilized to estimate it, as 
shown in Table 2. 

Si and S 2 at HA were significantly higher for both 
AVVAHR and V E /V C o 2 methods (Table 3). The ratio 
S 2 /Si tended to be higher at HA regardless of the 
estimation method utilized (Table 3), and this was par- 
ticularly so for the AV^/AHR method. 

Discussion and conclusions 

The main findings of this study are as follows: (1) VCP 
is clearly identifiable during incremental exercise at HA; 
(2) in conditions of persistent hypoxia, such as that elic- 
ited by HA exposure, VCP can be easily estimated by the 
AV^/AHR analysis. In the present study on young 
healthy individuals, applying standard methods for the 
identification, VCP could be identified both at SL and 
HA. Our data are in accordance with the observation of 
Agostoni et al, who demonstrated that at HA the VCP 
was clearly identifiable [22] . 

Of note, we observed a 30% reduction in Vo 2 -VCP 
when moving from SL to HA. The lower VCP at HA 



likely reflected the shortening of the isocapnic buffering 
period because of the increased hypoxic sensitivity of the 
carotid chemoreflex [22]. A further influence was pre- 
sumably the premature development of metabolic acid- 
osis during exercise at HA, as demonstrated by the 
lower 9 L likely due both to the reduction in arterial 0 2 
content and impaired convective 0 2 transport to the ex- 
ercising muscles [23]. 

Importantly, as demonstrated by the Bland-Altman 
analysis, we were able to demonstrate that the AV E / 
AHR approach is a reliable method for VCP estimation 
during incremental exercise at HA. The limits of agree- 
ment and the coefficient of variation between the Ve/ 
VcOi and AV^/AHR methods for VCP estimation are 
quite narrow and likely to fall within the between-day 
intra-subject Vo 2 variability [24]. 

The results of the present study are in agreement with 
a previous work of our group that demonstrated the reli- 
ability of the AV^/AHR breakpoint for VCP estimation 
at SL during incremental exercise [10,11]. At SL, assum- 
ing a linear relationship between Ve and HR (at least up 
to the heavy-intensity domain), the steepening of the 
AV^/AHR slope is more evident above VCP where (1) a 
steeper increase in Ve / 0 2 relationship is expected [9] 
and (2) the HR response could either maintain the same 
linearity as for more moderate exercise or in some 
instances (depending on exercise performance and 
protocol) reduce its rate of increment with respect to 
work rate [25]. Interestingly, at HA, we observed higher 
Si and S 2 AV^/AHR slopes. These findings are in ac- 
cordance with previous studies demonstrating that Ve at 
rest and at any given work rate during exercise is in- 
creased at HA because of the greater hypoxic drive 
[12-26] and also because of the reduction in the rate of 
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increase in HR between rest and peak exercise that is 
commonly observed at HA [15,16]. Moreover, compared 
to SL, we observed an increase in the S 2 /Si ratio for 
AVVAHR at HA. This finding suggests that S 2 is 
influenced both by hypoxic ventilatory drive and pos- 
sibly by an augmented contribution to V e from the 
metabolic acidosis above the VCP [1]. However, our rea- 
soning remains speculative as, for technical reasons, we 
were unable to measure arterial (or arterialized) [lactate] 
or pH during the exercise. 

Below VCP, an earlier (in the Si region) but less evi- 
dent breakpoint in AV^/AHR was discernible in some 
instances (5 out of 12), but only in one out five cases 
such a breakpoint was coincident with the anaerobic 
threshold. 

We are aware that our study has limitations, particu- 
larly with regard to the small sample size and the 
constrained characteristics of our study population, 
which prevented the evaluation of factors such as age, 
sex, and level of fitness, each of which is known to influ- 
ence the variables of interest. Also, we are aware (1) that 
the choice of a rapid work rate incremental protocol 
may have influenced the results utilized, i.e., different re- 
sults may be observed if a slower work rate increment 
protocol is used, and (2) of the lack of validation for pul- 
monary gas exchange measurement at HA. We took 
precautions to limit errors in the calibration procedures, 
with all equipment being calibrated before each test and 
all tests being performed indoors in the Pyramid Labora- 
tory at an ambient temperature of approximately 20°C 
and the necessary corrections for barometric pressure, 
humidity, and environmental temperature being applied 
through the dedicated software. Thus, our investigation 
should be considered as a pilot study conducted at ex- 
treme altitude, with the conclusions only being applic- 
able to healthy untrained young adults. 

In conclusion, we were able to demonstrate the reli- 
ability of the AV^/AHR method for the estimation of 
VCP at HA in a small group of sea-level residents, an 
environmental condition associated with different venti- 
latory and cardiac responses to exercise compared to SL. 
The attraction of the Al/g/AHR method is that it is a 
less expensive method for VCP estimation, compared to 
those that utilize expired gas measurements, and may 
therefore be utilized in sports medicine and in extreme 
conditions such as high altitude. 
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